Scalar chiral ground states of spin ladders with four-spin exchanges 
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We show that scalar chiral order can be induced by four-spin exchanges in the two-leg spin ladder, 
using the spin-chirality duality transformation and matrix-product ansatz. Scalar-chiral-ordered 
states are found to be exact ground states in a family of spin ladder models. In this scalar chiral 
phase, there is a finite energy gap above the doubly degenerate ground states and & Z2 'X Z2 x. Z2 
symmetry is fully broken. It is also shown that the SU(4)-symmetric model, which is self-dual under 
the duality transformation, is on a multicritical point surrounded by the staggered dimer phase, the 
staggered scalar chiral phase, and the gapless phase. 
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I. INTRODUCTION 

Recently, four-spin exchange interactions have been at- 
tracting interest in spin ladder models and spin-orbital 
models, because these interactions in fact appear in 
many system3ii2iiiiSiSiL& and can induce exotic ground 
stateSi^i-^^'-^^i-^^'-^'^i-^^'-^^i^^i-^'' Various types of four-spin in- 
teractions appear associated vifith diverse mechanisms, 
e.g., cychc-exchange processes Coulomb repulsion 
between doubly degenerate orbitals^ and spin-phonon 
couplings. Though two-leg antiferromagnetic Heisenberg 
spin ladders show a rung-singlet ground state and have an 
energy gap between the ground state and excitationspi 
it was revealed that four-spin exchanges can induce a 
gapped staggered dimer (or spin-Peierls) phasa2ii2iii and 
a gapless phaseii*i2iiiiiii^ around the SU(4)-symmetric 
point and that a gapped phase with a dominant vector 
chirality correlation also appearsiiSiii 

Very recently, Lauchli, Schmid, and Troyer^-^ numeri- 
cally found a new scalar chiral phase in the two-leg spin- 
1/2 ladder with four-spin cyclic exchange. Scalar chiral 
states, in which both the time-reversal and parity sym- 
metries are broken, had been discussed in the context 
of anyon superconductivity^^'^^ and the anomalous Hall 
eSectr^^ but realization of scalar chiral order in SU(2)- 
symmetric systems had been difficult and a challenging 
problem. A scalar chiral order due to the four-spin cyclic 
exchange was first proposed on the triangular lattice for 
magnetism of solid "^He filmsj^SiSi though finite-size sys- 
tem analysis could not find evidence for such ordering, 
instead showing spin-liquid ground statesi^SiS Study of 
spin ladders with four-spin exchanges is expected to clar- 
ify the possibility of exotic magnetism induced by the 
four-spin interactions. 

In this paper, we study the two- leg spin- 1/2 ladder 
with four-spin exchanges, whose Hamiltonian is given 
in the next subsection, and give a rigorous example of 
a scalar chiral ground state. In our analysis, the spin- 
chirality duality transformation we introduced in Ref. ^| 
plays an important role. A new class of models that have 



an exact ground state with scalar chiral order are con- 
structed by the means of matrix-product ansatz. The 
scalar chiral phase has the following nature: (i) the 
ground states are doubly degenerate, (ii) there is a finite 
energy gap between ground states and excited states, (iii) 
the ground states have long-range staggered scalar chi- 
ral order and exponentially decaying spin correlations, 
and (iv) a Z2 x Z2 x Z2 symmetry is fully broken. It is 
also found that the phase boundary of the scalar chiral 
phase touches the SU(4)-symmetric point and the SU(4)- 
symmetric model is on a multicritical point surrounded 
by the staggered dimer phase, the staggered scalar chiral 
phase, and the gapless phase. 

This paper is organized as follows. The model Hamil- 
tonian is given in the next subsection. In Sec. ^ we 
summarize the duality transformation and its applica- 
tion to the present model. In Sec. IIIII it is shown that 
a scalar chiral state is the exact ground state in a pa- 
rameter region of the model Hamiltonian. The nature 
of the scalar chiral phase is discussed. In Sec. IIVI we 
discuss the phase diagram around the SU(4)-symmetric 
point and find that the SU(4) model is on a multicritical 
point. Finally in Sec. we conclude with a discussion. 
Appendix contains a unitary description of the duality 
transformation. 



Model Hamiltonian 

The Hamiltonian of the two-leg spin-1/2 ladder with 
extended four-spin exchange interactions is defined as 

"H = J\^^isi,l ■ Si,i+i + S2,l ■ S2,l + l) + Ji'^Sij ■ S2.1 

I I 



Jd ^ (si,; ■ S2,;+i + S2,; ■ Si^;+i) 
•^rr ^ (si,; ■ S2,;) (si,;+i ■ S2,;+i) 
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+ -'"dd ^ {Sl.l ■ S2.I+1) (S2,i • Si,i + i) . (1) 
I 

This Hamiltonian includes a variety of models: (I) Four- 
spin cyclic-exchange model. When four-spin exchange 
constants satisfy Jn- = Jn = — Jdd, the four-spin terms 
describe the cyclic-exchange interaction.^ (II) SU(2) x 
SU(2) model. When the parameters satisfy — Jd = 
Jj.j. = J(jd ~ 0, the Hamiltonian has an SU(2)x SU(2) 
symmetry. This model was studied extensively as the 
SU(2) X SU(2) spin-orbital model. It was revealed 
that when Ju > the ground state has a staggered 
dimer (or spin-Peierls) orderiSiiS for AJ\ > J\i and is 
gaplessi2iiiiii^ for - Jn < 4Ji < Jn- (HI) SU(4) model. 
As a special case of model II, the Hamiltonian has an 
SU(4) symmetrjsiii2£ at AJ\ = J\\, which was exactly 
solved by the Bethe ansatzf^ 

II. DUALITY 

A. Duality transformation 

Let us begin with the spin-chirality duality transforma- 
tion, which we developed in Ref. We introduced the 
duality transformation defining new spin- 1/2 pseudospin 
operators 

S/ = ^ (Sl,i + S2,;) - Si^i X S2,i, (2) 
T/ = ^ (Si,i -f S2,;) -|-Si,i X S2,i, (3) 

which obey the commutation relations of spins and sat- 
isfy {Sf y — {Tj^Y = 1/4 for a — x,y or z. In the same 
way, the original spins Si^; and 82.; are expressed in terms 
of Si and T; in the forms Si,; ^ \ (S; + T,) -l-S; x T; and 
S2.; = i (S/ -I- T;) — S; X T;. In Appendix, we show that 
this transformation derives from a unitary operator J77r/2 

in the form S; = C/^/2Si, if/2/2 = f/^/2S2,zC/2/2- 

Since the following relations hold, 

Sl,i + S2,; = S; + Ti, 

sij -S2,i = 2 S, X Tz, 
-2 si,i X S2,i = Si - Ti, 

this transformation exchanges the Neel-type spin and 
vector chirality degrees of freedom on the same rung. As 
an example of the duality, one can show that the trans- 
formation of the dimer (or spin-Peierls) operator 

Od(0 = Si^i • Si,( + i - S2,i • S2,i + 1 (4) 

leads to the scalar chiral operator 

Csc(0 = (Sl,; + S2,i) • (Sl,i + 1 X S2,i + l) 

+ (si,i X S2,i) • (si^i+1 + S2,i + l). (5) 

It is known that a group of two-leg ladders with four- 
spin interactions shows the staggered dimer order in 



the ground stateiSiiSiii The above duality relation hence 
shows that their dual models have the staggered "scalar 
chiral" order in the ground state. We will discuss these 
models in Sees. IIIll and II VI 

We now consider the transformation of the spin states 
on rungs. Since the total spin on each rung is conserved 
under the transformation, the fully polarized spin state 
I T)i,i| 1)2, i is transformed to | t)s,;| T)t,/, where |a)^,i 
(a =t, I and [i = 1,2) denotes the spin state operated by 
s^ i, and \a)s.i {\ce)T.i) is the pseudospin state operated 
by Si (Ti), respectively. By applying Sf and Tf to 
I T)s,/| Dt,;, all pseudospin states in dual space can be 
constructed in the forms 

I T)s,/| T)t,; = I T)i,i| T)2,i, 

g-i7r/4 

|T)5,/U)T,i - ^^(1 T)l,iU)2,i +iU)l,i| T)2.i), 
iir/i 

I i)sA T)T,i - ^(1 T)i,n 1)2,1 -i\ T)2,0, 

I i)s,i\ Dta - I i)u\ i)2,l. (6) 

Using eigenstates for the total spin on each rung, one 
finds that this transformation corresponds to a gauge 
transformation of the singlet bond state \s)i ~* ~i\s)i, 
while it keeps the triplet states invariant. See also Ap- 
pendix for further arguments. 

B. Duality in the model Hamiltonian 

We now apply the transformation to the model and 
consider the duality relation in parameter space. Under 
this duality transformation, the total form of the Hamil- 
tonian Q remains invariant, but the couplings change. 
The couplings of the dual Hamiltonian Ti. are given by 





= Jr, 








Jl 


- ^(^'^ 


^ Jd)4 


liJ.- 


- Jdd): 


Jd 


- ^(^'^ 


-Jd)~ 


liJ.- 


- Jdd), 


Jn 


= 2(Ji- 


Jd) + 


liJu + 


Jdd), 


Jdd 


= -2(Ji 


-Jd) 


+ liJn 


+ Jdd) 



To see the mapping in the parameter space, we rewrite 
the Hamiltonian in the form 

W = Jr ^ Si^i • S2,i + Jrr ^ (Si,i • S2,i) (Si,i+i • S2,i+l) 
I I 

+ M/"^ (Slj + S2,i) ■ (Si,i+i -f S2,i + l) 
I 

+ ^y^{(si,i ■ Si^i + i) (S2,i ■ S2,i+l) 
I 

+ (Slj • S2J+1) (Si,i+i • S2,i)} 
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S2,i) ■ (Sl,( + 1 - S2,( + l) 



+4(si,i X S2,;) • (Si,; + i X S2^( + i)} 

^^{(Sl,i - S2,() • (Si^( + i - S2,i + l) 
/ 

-4(si,; X S2,i) • (Si,i + i X S2^( + l)}, 



(8) 



where 



W = ^(Ji + Jd), X=^(Ju + Jdd), 

Y = ^(Ju- Jdd) + J(Ji- Jd), 

Z = -^(Ju- Jdd) + 7(Ji- Jd). 
Id 4 



(9) 



Straightforward calculations show that the duality trans- 
formation maps the parameters {J^,J„,W^X^Y,Z) to 
( Jr, Jrr, W^, -'^j i^, — -Z^); the transformation changes only 
the coupling Z of the last term to — Z, but leaves the 
other terms invariant. The Hamiltonian |^ is, thus, self- 
dual at the surface defined by Z = 0, i.e.. 



Ju - Jdd - 4( Ji - Jd) = 0. 



(10) 



The parameter space of the Hamiltonian has six dimen- 
sions in total, and the self-dual surface divides the param- 
eter space into two regions Z > and Z < 0. It should 
be remarked that the SU(4)-symmetric model (4Ji = Jn, 
Jr = Jd = Jit = Jdd = 0) cxists On the sclf-dual surface. 
Duality around this specific model will be discussed in 
Sec.Hfl 



C. U{1) symmetry in the self-dual models 

Here we describe a U(l) symmetry in the self-dual 
models. Consider a continuous transformation with the 
following unitary operator 



Ue = J|exp 



This is a continuous extension of the duality transforma- 
tion (see Appendix) and it continuously transforms the 
dimer operator to the scalar chiral one. One can show 
that the Hamiltonian is invariant under this transforma- 
tion as UgTLUl = H for arbitrary 6 ii Z = 0. Thus the 
self-dual models are isotropic under the continuous rota- 
tion with the generator Si,; • S2,;, whereas the Z-term 
of the Hamiltonian JS)) lowers the symmetry. 



III. MODELS WITH EXACT SCALAR CHIRAL 
GROUND STATES 

In this section, we discuss an exact scalar chiral ground 
state of the Hamiltonian Q with the periodic boundary 



condition. To obtain the ground state 4'o, we use matrix- 
product (MP) states. We start from the following ansatz: 

l'o(u) = tr{gi{u)g2{-u) ■ ■ ■ g2N~i{u)g2Ni-u)}, (11) 

where m is a real variable and 

~ ' ^ 1 / iu\s)i + \to)i 
2 V %/2|t_i)i 



gi(u) 



-V2\t+i)i 
iu\s)i - \to)i 



(12) 



l-{iul\s)i - V2<j+\t+i)i + V2a-\t^i) 



\to)i}- 



Here \s)i and \t^)i are, respectively, the singlet and triplet 
states of the Ith rung, 2N is the total number of rungs, 
1 is the 2x2 unit matrix, and cr^ are the Pauli matrices. 
This form of the MP state can be obtained by the duality 
transformation of the MP state discussed by Kolezhuk 
and Mikeska,^" 

^km{u) = tr{gi(u)g2(-w) • • • 92W-i(")fl'2Jv(-'«)}, 

(13) 

where 

9i{u) = ^{ul\s)i - V2a+\t+i)i + V2a-\t^i)i + a,\to)i}. 

(14) 

It was shown that for several models this MP state 
^km(u) is the exact ground state with a staggered dimer 
order. At u = and u = oo, the two states ^^oiu) and 
^km(w) are equivalent, which means that each of ^'o(O) 
and 5'o(oo) is self-dual. 4'o(0) and \t'o(oo) are, respec- 
tively, an Affleck-Kennedy-Lieb-Tasaki (AKLT) state^S 
and a rung-singlet state. For < u < oo, however, ^o('^) 
and 5'km(m) are orthogonal in the limit N ^ oo. 

Since 4'km(w) (0 < m < oo) has the staggered dimer 
order, \E'o(m) has the staggered scalar chiral order because 
of the duality relation. Using the technique developed 
by Kliimper, Schadschneider, and Zittartz^S*^ one can 
evaluate the scalar chiral correlation in (u), 



{Oscil)Oscim)) - (-1 



12m 



(u2 + 3)2 



(15) 



for |Z — m| > 1 in the limit N oo and show spontaneous 
breakdown of the chiral symmetry 

12m 



(Osc(0> = (-1)' 



(m2 + 3) 



2 ■ 



One can also evaluate that ^'o(u) has no dimer correla- 
tion 

(Od(0Od(™)) =0. 

In the same way, the spin and vector chiral correlations 
are obtained as 



1 



- 1 



(u2 + 3)(u2-l) \u^ + 3 

((Sl,i X S2j)"(si,m X S2,m)") 



\l—m\ 



(m2 + 3)(m2-1) \^u2 + 3 



2\ 



(16) 



(17) 
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for a = X, y or z. The spin and vector chiral correlation 
lengths are equal and given by — = ln{(u^ + 
3)/|u^ — 1|}, whereas the scalar chiral correlation does 
not have any exponentially decaying term. 

By the duality transformation of the MP-solvable mod- 
els presented by Kolezhuk and MikeskajiS we find that 
the MP state 1)11(1 is an exact ground state of the fol- 
lowing three classes of models. One can prove that the 
state \l/o is a ground state in these models, reducing a 
local Hamiltonian on the Ith and I + 1th rungs 

to a positive semi-definite form — -Bo) ^ 0, where 

Ti. = and showing that ^'o has zero eigenenergy 

in the reduced Hamiltonian as (/if^j+i— i5o)fl';fl'/-i-i=0. The 
reader who is interested in the method of proofs should 
refer to Ref . E 

(A) Scalar chiral models. For a family of models 



Jd 
Jrr 



8J(2 - 3y) 4J(1 - y) 



3(4 -3y) ' 
J(8 - 9y) 



Jn = 



4-32/ 
16J 



= 



3(4-3?;)' 3(4 -3y)' 

-4J2; 



Jdd — 



3y' 



(18) 



with < y < 1 and J > 0, the ground states are dou- 
bly degenerate and given by ^'o(l) and \E'o(— 1). The 
ground-state energy per rung is Eq ~ —'iJ/'i. This 
model is dual to the "checkerboard-dimer model" given 
in Ref. 0, which has a staggered dimer order, and 
hence from the duality relation the present model be- 
longs to the scalar chiral phase. Excitations of the 
staggered dimer phase were studied by variational trial 
stateSfiS numerical calculations,^^ and field-theoretical 
analyses li^ii^ii^ These studies concluded that excitations 
have a finite energy gap. In fact, extending the argu- 
ments by Knabe;^^^ we can prove the finiteness of the 
energy gap in the checkerboard-dimer model dual to the 
model p8|l with y = 2/3 and in a finite region around 
this point. From the duality we conclude that, in the 
scalar chiral phase, there is a finite energy gap between 
the ground states and excited states. The analysis in 
the dual modeliS indicates that, sA y — 1, the present 
system enters the fully polarized ferromagnetic phase 
through a first-order transition. Furthermore, we can 
extend the parameter space which has the exact scalar 
chiral ground state. For y = 2/3, the dual Hamilto- 
nian has an SU(2) x SU(2) symmetry and the Hamil- 
tonian is written as a product of projection operators 
(4J/3)Ei(si,i • + 3/4)(s2^, • S2,i+i + 3/4). Then 

one can construct the model with the exact scalar chiral 
ground state by generalization of projection operators^I 
and the duality transformation. 

(B) Model at a phase boundary between the scalar chiral 
and staggered dimer phases. At y = of the model (|18|l , 



Ju 



4J 

T' 



Jl = J, 



2 J 

Jd = -J, Jrr = Jdd = 0, (19) 

with J > 0, the ground states are given by ^'o(w) with 
arbitrary u and highly degenerate. Note that this model 
is equivalent to the "multicritical model" in Ref. ^3 and 
self-dual under the duality transformation. The scalar 
chiral model 1)181) . thus, connects with the checkerboard- 
dimer model at this special parameter point. At this 
phase boundary, both one magnon and a pair of scat- 
tering solitons have energy gaps in the staggered dimer 
state, as discussed by Kolezhuk and Mikeska.^° However, 
because of the U(l) symmetry in this model, the gen- 
erator Si_i • S2,; can create gaplcss collective (Goldstone) 
modes, which are singlet bound states of two magnons. 
In fact, one can show that the following trial state be- 
comes gapless aX p — Q, IT {Q — —1) and p = 7r/2 (C = 1): 

(i-i 



l*o(p))L = ^e^^^'Tr m92-i(l)92.(-l)9^'^(l) 
X g2i+2{-l) n to-i(l)g2i(-l) > , 



i=l+2 



fT"92i-l(l)^7"92/(-l) W(l) 



C92;-i(l) 



92;+i(l) 



(C) Model with two second- order phase boundaries. For 
a family of models 



Jr 

Jl 

Ju 



-Jrr = ^(u2-l)(u2+3), 

6 



48 
J 

12 



(3u2 + 5)(m2 + 3), 
(Sw^-f 3)(m2 + 3), 



Jd = 77" 



Jdd = -t(w^ 
6 



(20) 
Qu^ - 3), 



with arbitrary u, the ground states are 5'o(u) and 
^'o(— u), and the ground-state energy per rung is Ea = 
-J{7u* + 22^2 + 19)/64. The model at u = 1 is equiv- 
alent to the model H18|l at y = 2/3. The arguments for 
the dual modeliS lead to the conclusion that the model 
1)21) |l undergoes a phase transition to the Haldane phase 
at u = and to the rung-singlet phase at u = cxd. Both 
of the transitions are of second order and accompanied 
with vanishing of energy gaps for solitons. 

In total, five phases appear in the MP-solvable models 
discussed above and in Ref. "T^. Some of the phase tran- 
sitions between them actually happen in the parameter 
space of the solvable models. These are summarized in 

Fig.m 

The nature of the scalar chiral phase is summarized as 
follows: (1) the ground states are doubly degenerate, (2) 
there is a finite energy gap between the ground states and 
excited states, and (3) the ground states have long-range 
staggered scalar chiral order and exponentially decaying 
spin correlations. 
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IV. AROUND THE SU(4)-SYMMETRIC POINT 
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FIG. 1: Schematic picture of five phases and phase transitions 
that appear in MP-solvable models. Phase transitions occur 
along arrows. The number attached to each arrow denotes 
the order of the phase transition. Black circle (MC) denotes 
the multicritical point 11911 . The scalar chiral phase is dual to 
the staggered dimer phase, while Haldane, rung-singlet, and 
ferromagnetic phases are self-dual. 



It is also easy to show that the string order, which was 
originally found in the Haldane state^S*^ exists in the 
scalar chiral state ^o(w). The expectation value is given 
by 

k-l 

+ 'h) n ^M^A^h + shmsi, + si,}) 

1=3 

= 4/(w2 + 3)2 (21) 

for a = x,y, or z. We note that the staggered dimer state 
^'km(u) also has exactly the same expectation value of 
the string order because the string operator is invariant 
under the duality transformation. This string order im- 
plies that a hidden Z2 x Z2 symmetry is spontaneously 
broken in the ground state4ii This hidden symmetry was 
foun d^9'i"l | by applying a nonlocal unitary transformation 
U, and one can find that this symmetry exists also in the 
Hamiltonian It should be noted that this Z2 x Z2 
symmetry is independent of the Z2 chiral symmetry as- 
sociated with the scalar chirality, since the scalar chi- 
ral operator after the nonlocal unitary transformation 
UOscU~^ has the corresponding Z2 x Z2 symmetry. In 
finite systems with open boundary conditions, the scalar 
chiral MP ground states in fact have eightfold degener- 
acy associated with boundary spins and chirality. Thus, 
the hidden Z2 x Z2 symmetry, as well as the Z2 chi- 
ral symmetry, is spontaneously broken in the scalar chi- 
ral phase. Recently, a different useful quantity Z2n — 
{exp[{2TTi/2N)Y^f^-^l{sl 1 + S2 i)]) was proposed,*^ which 
detects the average number n of valence bonds between 
neighboring rungs as limN^oo Z2N = (—I)"- In ^o(u)j 
the expectation value is estimated as limAr^oo Z2n = ^ 1 
for finite u. This is consistent with the above valence 
bond picture, because Z2n is also invariant under the du- 
ality transformation, and n = 1 in the staggered dimer 
state. 



Using the duality relation, we discuss the phase dia- 
gram around the SU(4)-symmetric point and show that 
this SU(4)-symmetric point is a multicritical point. 



A. SU(2) x SU(2) spin ladders 

We start from two SU(2) x SU(2) spin ladders. One 
SU(2) X SU(2) spin ladder is model II. Here we consider 
the case Ji > 0. For Ji/Jn = 1/4, this model is SU(4) 
symmetric and exactly solvable by the Bethe ansatz, and 
the ground state is gapless critical. ■^-^ For Ji/Jn > 1/4, 
the ground state has a finite gap and a staggered dimer 
Qj.(j^gj.^iOjiiji2 -^^j-^gj-gas for < J\/Ji\ < 1/4 the ground 
state is gapless and criticaliii^iSiiiiiii^ The total Hamil- 
tonian can be divided into the part of the SU(4) model 
'Ho (Jn — 4Ji with fixed Ji) and the perturbation V in 
the form 

n' = Ho + XV, (22) 

Ti-O = J\ y^{(sij ■ Sij+i -I- S2,; ■ S2,;+i) 
I 

-f4 (Sij • Sij+i) (S2J • S2J+1)}, 

V = ^ (Sl,i ■ Sl,i+l) (S2,i • S2J+1) , 

where A = — 4Ji -I- J\\ and J\ is fixed. Renormalization 
group analysis concluded that if the parameter A is nega- 
tive, the perturbation V is relevant and leads to a genera- 
tion of a staggered dimer order with a finite spin gap, and 
if the coupling parameter is positive, this perturbation is 
irrelevant and keeps the ground state gapless li^i^ii^ 

Applying the spin-chirality duality transformation to 
Eq. H22|) . one obtains the dual Hamiltonian 

ii' = Ho + XV, (23) 

V = g (si,; ■ Si^;+i -t- S2,; ■ S2,;+i) 

- o • S2,;+i -I- S2,; ■ Si^;+i) 

+ ^^{si,l-Si^l+i){s2J -82,1 + 1) 
I 

+ ^^{si,rS2,l + l){s2,l-Si,i + i). 



Here Ho is self-dual, and V is the perturbation dual to 
V. The couplings of H' , in total, are given as 



Jr = Jn=0, 



Jn = 2Ji + ijii, Jdd = -2Ji + ^ Jii. (24) 
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FIG. 2: Schematic possible pliase diagram around the SU(4)- 
symmetric point. The phase transition on the solid line is of 
second order and that on the dashed line is either first order 
or second order. Other phases might be inserted around the 
phase boundaries. 

From this transformation, one finds that the Hamiltonian 
H28|l has a hidden SU(2) x SU(2) symmetry, where gen- 
erators are given by J2i ^^'^ J2i '^t ct = y, or z. 
The duahty transformation leads to the case where, if the 
coupling parameter A is negative, the perturbation V is 
relevant and induces a staggered scalar chiral order with 
a finite spin gap, and if the coupling parameter is posi- 
tive, this perturbation is irrelevant and keeps the ground 
state gapless. When A — —\J\ (i-e., J\\ = ^J\)^ the model 
7i' equals to the scalar chiral model H18|) with y = 2/3 
and as we have shown in Sec. Ill the exact ground state 
has an energy gap and the scalar chiral order. 

B. Phase diagram around the SU(4) point 

We now discuss the phase diagram around the SU(4)- 
symmetric point, and the two SU(2) x SU(2) models 
given above, considering the following generalized Hamil- 
tonian: 

=Ho + AiV + A2V. (25) 

This Hamiltonian contains two kinds of perturbation to 
the SU(4)-symmetric model. Because of the duality, 
phase boundaries are symmetric with the line Ai = A2. 
The nature of phases in Ai > A2 is related to that in 
Ai < A2 by the duality transformation. The above con- 
sideration leads to a conclusion that the SU(4)-symmetric 
point is a multicritical point and surrounded by the stag- 
gered dimer phase, the staggered scalar chiral one, and 
the critical one. If the scalar chiral phase touches with the 
staggered dimer phase, the phase boundary between two 
phases must exist exactly on the self-dual line Ai = A2 
(see Fig. 121). Because of the U(l) symmetry, a rigorous 
theorem^ concludes that in general both orders disap- 



pear on the self-dual line and hence the transition be- 
tween the scalar chiral and dimer phases is second order, 
but, if uniform susceptibility of Si,; • 82,; diverges, both 
orders can exist on this line. Note that the latter actually 
happens in the model (|19|l . One plausible phase diagram 
is shown in Fig. 13 Recently we have numerically studied 
this phase diagram and obtained results consistent with 
the present conclusioni^ 



V. DISCUSSION 

In this paper we have shown a rigorous example of 
scalar chiral ground states in SU(2) spin ladders with 
four-spin exchanges. The exact duality relation is the 
keystone of our theory. Our results demonstrated that 
four-spin exchanges can actually induce the scalar chi- 
ral long-range order. The scalar chiral phase extends to 
a wide parameter region and touches with the SU(4)- 
symmetric point. Previously, a scalar chiral phase was 
numerically found in the four-spin cyclic exchange model 
on the two-leg ladder.^^ In their phase diagram, the scalar 
chiral phase appears next to the staggered dimer phase 
and the phase boundary indeed exists on the self-dual 
point J 4^1 J = 1/2. This situation in the vicinity of the 
self-dual point shows a resemblance to that around the 
self-dual line (Ai = A2 < 0) in Fig. [3 Our recent numeri- 
cal study of the Hamiltonian Q indicates that the scalar 
chiral phase we found in this paper extends to the four- 
spin cyclic-exchange case and that two phases belong to 
the same one. 

We have shown that the SU(4)-symmetric model is 
self-dual under the spin-chirality duality transformation. 
We here note that this statement holds for the SU(4) 
spin-orbital models on arbitrary lattices. Recently SU(4) 
spin-orbital models on two-dimensional lattices^ and on 
laddergi^ have been studied and it was discussed that a 
plaquette ordering may appear in the ground state. On a 
four-site plaquette, the SU(4) singlet state is the unique 
ground state and therefore it must be self-dual under the 
duality transformation. We hence conclude that plaque- 
tte ordering is also self-dual. 

Last, we discuss the universality classes of phase tran- 
sitions. The phase transitions into the scalar chiral phase 
are naturally in the same universality class as the dual 
transitions into the staggered dimer phase. For example, 
since the phase transition between the rung-singlet phase 
and the staggered dimer phase belongs to the c = 3/2 
SU(2)2 criticality,- we conclude that the transition be- 
tween the scalar chiral phase and the rung-singlet phase 
also belongs to the same one. Since the two-dimensional 
Ising model is related to the c = 1/2 criticality, this 
c = 3/2 criticality can be plausibly regarded as a con- 
sequence of the ^2 X ^2 X Z2 symmetry breaking. 
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APPENDIX: UNITARY OPERATOR FOR 
DUALITY TRANSFORMATION 

In this appendix, we show that the duality transforma- 
tion ||2J) and (PI corresponds to a unitary transformation 
of two spins on rungs. The unitary operator is given by 



Ug = Y[exp[-iePi{s)] 



J|exp 



(26) 



with = 7r/2, where Pi{s) denotes the projection opera- 
tor onto the singlet state on the lih rung. Since the gener- 
ator is transformed as Si; -82^; — (si_i-t-S2,i)^/2— 3/4, this 
unitary conserves the total spin on each rung. Note that 
the generator of this unitary is a summation of SU(4) gen- 
erators s"s2 ("^ ~ ■^i y^ hence the SU(4)-symmetric 
model is naturally invariant under this transformation. 

Let us demonstrate the unitary transformation of 
spins. It is convenient to reduce the unitary operator 
to the form 



(27) 



where we have used the relation [P;(s)]^ = Pi{s). Using 
the following commutation relations 

[si,i ■ S2,i, si^i + S2.1] = 0, (28) 
[si,i ■ S2,i,si^i - S2.1] = 2isi^iXS2.i, (29) 

1 

[Sl,l ■ S2,l,Sij X S2,l] = ~-^isi,i - S2,l), (30) 

and the unitary relation UgUe^ — 1, one can perform the 
unitary transformation of spins in the forms 



(31) 



Ueisij + S2j)Ub^ = si^i + S2,i, 

Ug{si^l - S2,l)Ue'' 

= cos6{si,i - S2,i) - 2sin0{si^i x S2,i), (32) 
Ueisij X S2j)Ue'^ 

= i sin6'(si,i - S2d) + cos9{si^i x S2,i). (33) 
When 6 — tt/2, we obtain the original duality transfor- 
mation 



Si — U^I2Si^iUt,/2\ 



(34) 
(35) 



From the form of the unitary operator, it is clear that 
this unitary corresponds to a gauge transformation of 
the singlet bond state 

Ub\s)i = e-''\s)u (36) 
Ue\t„^)i = |t™);, (m = -1,0,1). (37) 

When 9 — 7r/2, one obtains the relation ^ from 

W)s,iW)t,1 = U^/2W)lM')2,l (38) 

for a (a') =t,i. 
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